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ABSTRACT 

We present the discovery of 5 mihisecond pulsars found in the mid-Galactic latitude 
portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have 
rotational periods from ^ 2.3 ms to ^ 7.5 ms, and all are in binary systems with 
orbital periods ranging from ^ 0.3 to ~ 150 d. In four of these systems, the most 
likely companion is a white dwarf, with minimum masses of ~ 0.2 Mq. The other 
pulsar, J1731— 1847, has a very low mass companion and exhibits eclipses, and is 
thus a member of the "black widow" class of pulsar binaries. These eclipses have 
been observed in bands centred near frequencies of 700, 1400 and 3000 MHz, from 
which measurements have been made of the electron density in the eclipse region. 
These measurements have been used to examine some possible eclipse mechanisms. 
The eclipse and other properties of this source are used to perform a comparison with 
the other known eclipsing and "black widow" pulsars. 

These new discoveries occupy a short-period and high-dispersion measure (DM) 
region of parameter space, which we demonstrate is a direct consequence of the high 
time and frequency resolution of the HTRU survey. The large implied distances to 
our new discoveries makes observation of their companions unlikely with both current 
optical telescopes and the Fermi Gamma-ray Space Telescope. The extremely circular 
orbits make any advance of periastron measurements highly unlikely. No relativistic 
Shapiro delays are obvious in any of the systems, although the low flux densities would 
make their detection difficult unless the orbits were fortuitously edge-on. 

Key words: pulsars: general - stars: neutron - pulsars: individual: PSR J1125-5825 - 
pulsars: individual: PSR J1708-3506 - pulsars: individual: PSR J1731-1847 - pulsars: 
individual: PSR J1801-3210 - pulsars: individual: PSR J1811-2405 



1 INTRODUCTION 

Millisecond pulsars (MSPs) are neutron stars (NSs) with 
rapid rotation rates that are believed to be formed in bi- 
nary systems when the NS accretes matter from the com- 
panion, causing a 'spin-up' in the NS's rotation rate (e.g. 
Alpar et al.|1982[ ). This process results in pulsars with spin 
periods between 1 and 100 ms and magnetic field strengths 
less than 10^" G. Those systems with spin periods less than 
about 20 ms are thought to have had low mass companions, 
while those with periods between 20 and 100 ms are thought 
to have been spun-up by a heavy white dwarf (WD) or NS 



companion. During the spin-up phase, the heating caused by 



accretion leads to the emission of X-ray radiation (David- 
[son fc Ostriker||1973[ ). These accreting systems are known 
as high-mass and low-mass X-ray binaries depending upon 
the companion mass (HMXBs and LMXBs respectively, see 
^Bhattacharya fc van den HeuveT] ( |1991[ ) for details of evo- 
lution). The link between LMXBs and MSPs was not con- 
firmed until pulsations with a period of 2.4 ms were observed 
in the accreting X-ray binary, SAX J1808.4— 3658 ( [Wijnands 
& van der Klis|1998| . This link has more recently been rein- 



forced by the discovery of PSR J1023-^0038 ( |Archibald et al. 
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20091, where radio emission is likely to have only switched 



on recently after an LMXB phase. 

In contrast to this evolutionary model, however, around 
20% of the MSP population appear to be isolated bodies, 
including the first M SP to be discovered, PSR B1937+21 
(Backer et al. 19821. Therefore, if this scenario is correct. 



rotating known pulsar is PSR J1748— 2446ad (Hessels et al. 
2006), which rotates at 716 Hz; the discovery of an MSP 
rotating even faster than this could rule out some NS equa- 
tion of state models (Lattimer & Prakash 2007 20101. By 



constraining the maximum mass of a NS, it should also be 
possible to eliminate many models of the equation of state 



and the isolated MSPs are the descendants of the eclipsing 
MSPs, we might expect some of them to eventually ablate 
their companions until nothing remains ( jRuderman et al.| 
1989). The discovery of the so-called 'black widow' pulsar, 
PSR B1957-I-20 (i Fruchter et al.|1988|) and other si milar sys- 
tems (e.g. PSR J2051-0827, [Stappers et al.|1996| , provides 
some evidence of this process taking place. 

In these systems, the pulsar's companion is typically of 
very low mass (~ 0.02 M©), and the pulsar is eclipsed by 
the companion for at least 10% of the orbit. As the pul- 
sar approaches eclipse, the pulses experience a delay due to 
the additional ionised gas through which the radiation has 
to pass, while outflowing material from the companion often 
gives rise to anomalous eclipses away from superior conjunc- 
tion. However, the timescale over which the ablation process 
would lead to the companion being destroyed is far too long 



(Lattimer & Prakash 2004'). To do this, one can use binary 



effects such as Shapiro delay which allow the pulsar mass to 
be measured (e.g. Lync et al.^2004i ); an effect which is most 
measurable for those systems where the pulsar is found to 



have a massive companion or a near edge-on orbit (e.g. De- 



morest et al. 2010 1. Therefore, the discovery of more systems 
where the NS mass can be measured (e.g. Freire et al.|2010" l 



would contribute to the understanding of the NS equation of 
state. For some pulsars in binary systems, it has been pos- 



sible to observe the companion optically (e.g. Bassa et al. 
2006 ) and measure its radial velocity as a function of orbital 
phase. However, for WD companions, often the faintness of 
the companion only allows the measurement of a tempera- 
ture, from which a mass must be estimated using evolution- 
ary models ( van Kerkwijk|[T996 1. 

Currently, pulsar timing arrays ( Hobbs et al.|20d9 Fer 



toexplain the abundance of isolated MSPs ( [Stappers et al.| |dman et al.|20l'ol |Jenet et al.|2009t are attempting to detect 



1996|). Clearly other mechanisms such as tidal disruption are gravitational waves by the correlation of arrival times (e.g 



needed to explain the existence of this group 



The discovery of PSR J 1903-1-0327 (Champion et al. 



penet et al. 2006 1 from pulsars. These arrays require long- 



2008 1 further challenges the conventional formation scenar- 



ios. Not only is the pulsar's companion a main-sequence star, 
but the orbit is also highly eccentric (e = 0.44), neither of 
which are predicted by the standard formation scenarios. 
Three alternative theories for the formation have been pro- 
posed; a) the pulsar was born in an eccentric orbit, spin- 
ning rapidly; h) the pulsar was spun-up in a globular cluster 
before being ejected into the Galactic disk; and c) the recy- 
cling of the pulsar in a triple system ( Champion et al.|2008 1. 
[Freire et al.| ( |2010[ ) have since confirmed that the main se- 
quence star is the companion, and excluded the possibility 
that PSR J1903-I-0327 is currently a member of a triple sys- 
tem, but maintain that the pulsar was born in a triple sys- 
tem, from which the component originally responsible for 
the spun-up of PSR J 1903-1-0327 has been ejected. However, 
if the alternatives are viable, we need to understand how 
many MSPs form via this channel and to do that a larger 
sample of MSPs is required. 

Beyond improving our understanding of binary star evo- 
lution, there are many areas of pulsar science which ben- 
efit from the discovery of recycled pulsars. The discovery 
of binaries with, for example, a neutron star or black hole 
companion offer the chance to test General Relativity (GR) 
and other theories of gravity in the strong field regime (e.g. 
Kramer et al.||2004l|2006||. The case of a pulsar with a WD 



term, high precision timing observations with high signal-to- 
noise ratio; MSPs, with their short spin periods and stable 



rotation rates, are well-suited for these timing arrays (Ver 



biest et al. 20091. However, not all MSPs have the narrow 



pulse profiles, regular rotation and high fiux density required 
to obtain high timing precision, nor is the distribution of 
the known sources as uniform over the sky as desired (see, 
e.g. [Hellings fc Downs] 1983[ ). Hence, further discoveries from 
pulsar surveys may make important contributions to pulsar 
timing arrays and improve their sensitivity to the stochastic 
background of gravitational waves. 

The largest previous pulsar survey, that found many 
of the known MSPs, was the Parkes Multibeam pulsar sur- 
vey (PMPS) by [Manchester et al.] ( |2001[ ), which surveyed 
the area of the Galactic plane bounded by 260° ^ / 5S 
50°,|fe| ^ 5°, discovering 26 MSPs. At higher Galactic lat- 
itudes, the Swinburne Intermediate-latitude pulsar survey 
(260° s; I s; 50°, 5° s; |6| ^ 15°), using a similar observ- 
ing system except for a reduced sampling time, discovered 
a further 8 MSPs (Edwards et al. ^OOlj [Edward s fc Bailes| 



2001 1. However, these surveys were limited by the relatively 



companion also offers the possibility of observing effects pre- 
dicted by GR. For example, PSR J1909-3744 is in a 1.5 day 
orbit with a companion of mass 0.2 M©, where Shapiro de- 
lay has been observed, allowing the measurement of both the 
companion and pulsar masses ( [Jacoby et al. [[2003] [2005] ). 

The NS equation of state is currently poorly under- 
stood, and the predictions of NS masses and radii vary dra- 



matically for different models (Lattimer & Prakash 20041 



One way to place a limit on these parameters is to iden- 
tify the limiting rotation frequency, beyond which the neu- 
tron star would break apart. Currently, the most rapidly- 



broad frequency channels (3 MHz) and coarse time sampling 
(250 /is and 125 /is, respectively) compared with today's dig- 
ital backends, which provide increased sensitivity to distant 
and rapidly-rotating pulsars. 

The ionised interstellar medium disperses the pulses 
emitted by radio pulsars as they traverse it. When remov- 
ing this dispersion, increased frequency resolution allows this 
correction to be made with reduced smearing of the pulses. 
The increased time resolution provides greater sensitivity 
to short-period pulses, such as those from MSPs, and by 
sampling with 2 bits (compared to 1 bit per sample in the 
PMPS) sensitivity is further increased by reducing losses due 
to digitisation. The High Time Resolution Universe Survey 
(HTRU, [Keith et al.[[20To| ) aims to make use of these im- 
provements in backend technology to perform a survey of 
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Table 1. Observational parameters for the mid-latitude portion 
of the HTRU survey. 



Number of beams 
Polarizations /beam 
Centre Frequency 
Frequency channels 



13 
2 

1352 MHz 
1024 X 390.625 kHz* 



Galactic longitude range 
Galactic latitude range 
Sampling interval 
Bits/sample 

Observation time/pointing 



-120° to 30° 
\b\ ^ 15° 
64 ^is 
2 

540 s 



*154 of these channels are then masked to remove interference 

the entire southern sky with much-improved sensitivity to 
MSPs. 



The ongoing PALFA survey at Arecibo (Cordes et al 



2006[ ) and the GBT 350-MHz survey ( |Boyles et al.||2010 l 

have each discovered several MSPs which occupy the short- 
period and high-DM region of parameter space which we 
look to probe with HTRU. In particular, the discovery, as 
previously mentioned, of PSR J1903+0327 ( [Champion et al.] 
2008 1 as part of PALFA stands out with a rotation period 
of 2.15 ms and a DM of 297.5 cm~^ pc as evidence that time 
resolution of 64 fis allows this parameter space to be probed. 

In this paper we outline the spin and orbital parameters 
of five MSPs discovered as part of the HTRU survey, com- 
pare the properties of these pulsars to the previously-known 
population, and study in detail the eclipses that one of them 
displays. These discoveries have been made with ~ 30% of 
the survey region observed, indicating that we might expect 
tens of MSPs to be discovered when the survey is completed. 



2 OBSERVATIONS 

2.1 Discovery and Timing 

The HTRU survey | |Keith et ar]|2010| ) is broken into three 
distinct components with different scientific objectives. The 
low latitude region covering Galactic latitude |fe| < 3.5°, the 
mid-latitude region covering |fe| < 15°, and the high lati- 
tude region which covers the remaining sky below declina- 
tion -1-10°. The mid-latitude survey is designed to find mil- 
lisecond pulsars that were missed by previous surveys due to 
excessive dispersion smearing but are still bright enough to 
be found in ~ 10 min. Observations are made using the 13- 
beam multibeam receiver on the 64-metre Parkes radio tele- 
scope, which has a half-power beam width of ~ 0.23°. Data 
were then processed using the hitrun pipeline described 
in [Keith et al.| all five MSPs were discovered in the mid- 
latitude component of this survey, which covers an area of 
the Galactic plane bounded by -120° < I < 30°,|6| ^ 15° 
with observations 540 s in duration. Full details of the survey 
parameters are given in Table [T] 

After their discovery and confirmation, timing obser- 
vations of these MSPs were made with baselines ranging 
from 250 days to more than 400 days. Those which are visi- 
ble from Jodrell Bank Observatory — PSRs J1731-1847, 
J1801-3210 and J1811-2405 — were regularly observed 
with the 76-metre Lovell Telescope in a band centred at 
1524 MHz, while PSRs J1125-5825 and J1708-3506 were 



Table 2. Observing system details for the timing observations 
made as part of this work. Note the specifications for the Lovell 
Telescope take into account the removal, as standard, of a section 
of the observing bandwidth. 



Telescope 


Centre Freq. 


BW 


■^chans 


(tobs) 




(MHz) 


(MHz) 




(s) 


Parkes 64-metre 


732 


64 


512 


900 




1369 


256 


1024 


600 




3094 


1024 


1024 


900 


Lovell Telescope 


1524 


384 


768 


900 


WSRT 


345.625 


80* 


512 


2400 



* With coherent dedispersion 

observed with the 64-metre Parkes radio telescope in a band 
centred at 1369 MHz (these Jodrell Bank and Parkes data 
are, hereafter, 'the 20 cm band'). At Jodrell Bank, observa- 
tions were made approximately once per week, and Parkes 
observations were more sporadic, with gaps between obser- 
vations sometimes lasting as long as three months. Timing 
observations at both observatories were made using digi- 
tal filterbanks (DFBs). Flux-calibrated observations using 
the Parkes DFB were used for measurement of the source 
fluxes. Calibration was performed using routines included in 
PSRCHIVE ( Hotan et al.|2004 | and with observations of both 
a calibrator source (Hydra A) and a pulsed cal. 

In order to study the pulse profile evolution with fre- 
quency, observations were made of some of the sources us- 
ing the 10/50 cm band receiver mounted on the Parkes 
telescope (the bands have centre frequencies of 732 and 
3094 MHz), and using the Westerbork Synthesis Radio Tele- 
scope (WSRT) and the PuMA II backend (Karuppusamy 
et al.|p008 l in a band centred at 347 MHz (hereafter, 'the 
92 cm band'). Observations were made with the WSRT for 
each of the pulsars with 5 > —35°, however, neither PSR 
J1801-3210 nor PSR J1731-1847 have been detected at 
this wavelength. The system parameters for all the observa- 
tions are shown in Table [21 

The timing solutions for each of the systems are shown 
in Table [3] The first quantities given in this table are the 
positions of the sources after timing over the specified range 
of dates. This timing has allowed the positions to be fitted 
to sub-arcsecond precision, with positional errors as given. 
Four of these new discoveries lie within 5° of the Galactic 
plane, with only PSR J1731 — 1847 outside this range, with 
& = 8.15°. 



2.2 Orbital and Spin Parameters 

The pulse periods of the newly discovered pulsars range from 

2.3 ms (PSR J1731-1847) to 7.5 ms (PSR J1801-3210) (see 
Table |3|. The change in period seen in confirmation observa- 
tions of all five sources made it apparent that the pulsars are 
in binary systems. With subsequent observations, we were 
able to fit for the observed variation in pulse period with 
time and deduce a set of orbital parameters for each pulsar. 
In order to precisely measure the spin and orbital parameters 
of each pulsar, the regular timing observations in the 20 cm 
band were used to produce a set of time-of-arrival (TOA) 
measurements. Using the standard pulsar timing procedure 
(e.g. Lorimer & Kramer 20051, these TOAs were used to 
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Table 3. Observed and derived parameters for the new millisecond pulsars. The DM distance has been estimated using the model of 
[Cordes fc Lazio] (|2002|l, while a pulsar mass of 1.4 Mq has been assumed in calculating companion masses. 



Parameter 


J1125— 5825 


J1708-3506 


J1731-1847 


J1801— 3210 


J1811— 2405 


Right Ascension (J2000) 


11:25:44.3654(3) 


17:08:17.623(1) 


17:31:17.6072(2) 


18:01:25.8890(2) 


18:11:19.8539(2) 


Decimation (JzUUU) 


cro OCT AC C7 f ^\ 

— 5»:25:1d.»d7(5) 


— 35:06:22.85(7) 


— 18:47:32.74(1) 


— 32:10:53.72(1) 


-24:05:18.72(7) 


Galactic Longitude (°) 


291.89 


350.47 


6.89 


358.92 


7.07 


Galactic Latitude (°) 


2.60 


3.12 


8.15 


-4.58 


-2.56 


T~\: o /1\T 

Discovery b/lN 


11 


18 


14 


11 


13 


Offset from survey 


0.12 


0.05 


0.06 


0.03 


0.17 


beam centre (°) 












TOA Range (MJD) 


55135-55431 


55127-55461 


55148-55399 


54996-55409 


55131-55389 




o.lUzzloyi6yoU4(f5 J 


4.0Uolooy4ozo(D ) 


Z.o4400y0ooo / / (Dj 




z.DDUoyoolDyU(l } 


P (xlO-20) 


5.963(5) 


2.3(4) 


2.49(1) 


0.265(97) 


1.36(2) 


DM (cm ^ pc) 


124.81(5) 


146.8(2) 


106.56(6) 


176.7(4) 


60.64(6) 


DM Distance (kpc) 


2.6 


2.8 


2.5 


4.0 


1.8 


Si. 4 GHz (mJy) 


0.44(2) 


0.33(1) 


0.60(2) 


0.21(1) 


0.37(1) 


Ll.4GHz (mJy kpc^) 


3.0 


2.6 


3.8 


3.4 


1.2 


T-c (yr) 


8.2 X 10* 


3.1 X 10^ 


1.5 X 10^ 


4.5 X 10^0 


3.1 X 10^ 


-Bsurf (G) 


4.3 X 10* 


3.2 X 10* 


2.3 X 10* 


1.4 X 10* 


1.9 X 10* 


i>lc (,G) 


1.3 X 10"^ 


3.2 X 10^ 


1. * X iU 


3.1 X 10^ 


9.3 X 10'* 


E (erg s"-*-) 


7.9 X 10^* 


9.9 X 10^3 


7.6 X 1034 


2.5 X 10^2 


2.8 X 10^4 


E/<i? (erg kpc-2 s'^) 


1.2 X 10^* 


1.2 X 10^3 


1.2 X 10^4 


1.6 X 10^1 


8.6 X 10^3 


Binary Model 


ELLl 


ELLl 


ELLl 


ELLl 


ELLl 


Orbital Period (d) 


76.4032169(4) 


149.13318(1) 


0.311134124(4) 


20.7716995(3) 


6.27230204(2) 


asini (It-s) 


33.638356(1) 


33.584233(5) 


0.1201594(9) 


7.809320(5) 


5.7056623(8) 


TASC (MJD) 


55126.348808(1) 


55132.23096(2) 


55132.3113000(8) 


55001.934481(3) 


55136.1686237(1) 


ei 


-0.00025349(5) 


0.0000004(2) 


0.00010(1) 


-0.000003(1) 


0.0000017(4) 


£2 


-0.00004404(6) 


-0.0002445(2) 


0.00000(1) 


0.000000(1) 


0.0000011(4) 


e 


0.00025720(5) 


0.0002445(2) 


0.00014(1) 


0.000004(1) 


0.0000024(4) 


u] (deg) 


260.140(2) 


179.91(5) 


264(5) 


290(20) 


253(7) 


Min. TTic (M0) 


0.26 


0.16 


0.037 


0.15 


0.24 


Med. mc (M©) 


0.30 


0.19 


0.043 


0.17 


0.28 


RMS of fit (/xs) 


5.538 


7.771 


6.339 


29.050 


5.539 


of fit 


0.95 


1.3 


7.8 


1.3 


1.3 



fit for each of the parameters presented in Table [3| using et al. 19991; and if the Shklovskii contribution was as 



the Tempo2 pulsar timing package ( Hobbs et al.|20d6| . The 
locations of these pulsars in the P-P plane are shown in 
Fig. [ij and are seen to lie towards the lower end of the pe- 
riod distribution, with each of the new discoveries fitting in 
the region traced by the previously-known population with 
short spin periods. 

The transverse velocity, Vr, of a pulsar at distance d 
makes a contribution to the measured value of P in what is 
known as the Shklovskii effect (Shklovskii jl970[ ) 

P IV^ 



P c d 

We can then calculate the value of Vt, given some con- 
tribution to P/P. Assuming this contribution is 10% for 
each MSP, we find that PSRs J1708-3506 and J1731-1847 
would be required to have very large values of Vt, 670 and 
350 kms~^ respectively. The other three MSPs would re- 
quire Vt < 100kms"\ which are 'realistic' values (^Toscanol 



high as 50% for these three, Vr would only need be ~ 
100 km s"^ It seems feasible, therefore, that there is a sig- 
nificant Shklovskii contribution to the measured value of P 
for these pulsars, which in turn implies that the values of 
E cc P and Bsurf oc P^/^ are overestimated and Tc oc P~^ 
underestimated in Tabled 

In this work we have used the ELLl pulsar timing model 
(as outlined in the appendix of [Lange et al.||2001[ ) to fit 
for all the relevant spin and binary parameters in each of 
the pulsars. This model is suitable for modelling the orbits 
of low-eccentricity systems, and for PSRs Jl 125— 5825 and 
J1708— 3506, gives a statistically significant value for the 
orbital eccentricity. With the exception of PSR J1731 — 1847, 
the reduced ~ 1 for each of the timing models. In the case 
of PSR J1731-1847, the higher can be attributed to the 
extra delay induced by the eclipses (as discussed below). 
With further observations, and taking care to remove TOAs 
around the eclipse region, the value is expected to reduce. 
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Binary Pulsars X 

Eclipsing Binaries 9 
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0,01 
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Figure 1. P-P diagram of the millisecond pulsars, divided by bi- 
nary 'type'. Lines of constant E cross the population at (from left 
to right) 10^^, 10^'^ and lO^'^ergs"^. Pulsars in globular clusters 
have been excluded from this figure. 



The orbital parameters for the five MSPs are shown in 
Table [3] from which it can be seen that they have a wide 
range of orbital periods, from 0.3 d (PSR J1731-1847) to 
150 d (PSR 1708-3506). Using the orbital parameters for 
each pulsar, we can relate the mass of the pulsar, nip, and 
the mass of its companion, rric, to the projected semi-major 
axis of the orbit, x, using the Keplerian mass function 



/(mp,mc) = 



(mcsini)^ _ in^ 



(2) 



where G is Newton's gravitational constant and Pb is the 
orbital period of the pulsar in the binary system, and i the 
inclination of the orbit. Assuming a pulsar mass of 1.4 Mq, 
the mass function may be solved for the mass of the compan- 
ion as a function of orbital inclination. Solving for i = 90° 
gives the minimum mass of the companion, and i — 60° 
gives the median mass; each of these values are given in 
Table ID Apart from PSR J1731-1847, the companions to 
these pulsars have median mass values of ~ 0.2 Mq, indicat- 
ing that they are typical pulsar- WD binary systems. PSR 
J1731— 1847, however, has a companion with a very low me- 
dian mass, only 0.043 Mq. As this low mass suggests (e.g. 



Ryba fc Taylor 1991 Stappers et al. 1996 1, PSR J1731-1847 



displays eclipses, which are discussed in detail in section |4] 



2.3 Pulse Profiles 

The pulse profiles of each pulsar in the different observ- 
ing bands are shown in Fig. [2] PSRs J1125— 5825 and 
J 181 1—2405 each display an interpulse, trailing the main 
pulse by ~ 0.55 and ~ 0.45 in pulse phase, respectively. 
The spectral properties of the two interpulses is different at 
shorter wavelengths — the interpulse of PSR J1125— 5825 
is not visible in the 10 cm band, probably due to the low 
S/N of the profile. However, in PSR J1811-2405 the flux 
of the interpulse increases relative to the main pulse at 
10 cm, indicating that the interpulse has a flatter spectrum 
than the main pulse. This has been observed previously, 
both for normal pulsars (e.g. Biggs et al.||1988 \ and MSPs 
J2322 -H2057 ( [Nice et al.|[l993l), J1012-H5 307 ( [Kramer et al.| 
1999[ ) and B1855+09 ( |Kijak et al!]|1997[ ). The interpulse of 



PSR J 1125— 5825 also appears to be complex, consisting of 
at least 2 components in the 20 cm band, while the profile 
of PSR J1811— 2405 displays ofT-pulse emission trailing the 
interpulse that can be seen in the pulse profiles at 10 and 
20 cm. 

The main peaks of the pulse profiles of PSRs 
J1731— 1847 and J1811— 2405 are narrow in comparison to 
the other new discoveries. This property can make high- 
precision timing of an MSP an easier task, however the 
eclipses of PSR J1731— 1847 make it unlikely to be suitable 
for pulsar timing arrays. The profiles of PSRs J1708— 3506, 
at 20 cm, J1801-3210, at 50 cm, and J1811-2405, at 92 cm, 
display broadening of the pulse with respect to the profiles 
at shorter wavelengths. In these cases, this broadening can 
be attributed to the scattering of pulses by the interstel- 
lar medium, as the high frequency resolution of the observ- 
ing systems precludes the possibility that this is caused by 
residual DM smearing across a frequency channel, while the 
observation in the 92 cm band has been coherently dedis- 
persed. 

Using the empirical fit of Bhat et al. (20041, we would 



expect the scattering tail in these observations to have a 
1/e timescale of 0.7 ms for PSR J1801-3210 and 0.1 ms for 
PSR J1811-2405, but only 0.02 ms for PSR J1708-3506. 
Since there is a large amount of scatter around the Bhat et 
al. relationship, the fact that this timescale is slightly under- 
estimated in the first two cases each case is not likely to be 
significant — for comparison, the NE2001 electron distribu- 
tion model ( Cordes fc Lazio|2002 I predicts scatter broaden- 
ing times of 0.01, 0.4 and 0.5 ms for PSRs J1708-3506, 
J1801-3210 and J1811-2405 respectively It is also pos- 
sible that this pulse broadening is due to intrinsic evolu- 
tion of the pulse profile with observing frequency, and for 
PSR J1708— 3506, it is unclear whether this broadening is 
due to a second component which trails the main peak. This 
component may be present, though faint, in the 10 cm ob- 
servation. 

For the two pulsars which were observed, but not de- 
tected, at 92 cm, assuming the standard relationship for 
pulse scattering as a function of wavelength (e.g. Lambert 
fc Rickett|2000| , oc A'"/^', the 1/e timescale for scatter- 
ing, Ts, should be increased by a factor of ~ 11 at 92 cm 
compared to 50 cm. In the case of PSR J1801— 3210, it 
seems likely that the increased scattering could be respon- 
sible for smearing the pulse over more than a pulse period. 
For PSR J1731— 1847, however, the observations were made 
well away from orbital phase 0.25 (removing the possibility 
of eclipse), nor does it appear that the scattering would be 
severe enough to smear the pulse by more than a pulse pe- 
riod, making it likely that the flux density at this wavelength 
was below the detection threshold, implying a relatively fiat 
spectral index, or a spectral break. 



2.4 Potential measurement of relativistic effects 

In certain cases, with favourable orbital parameters, it is 
possible to measure the mass of both the pulsar and its com- 
panion (e.g. [Demorest et al.||2010[ ). When light passes near 
a massive body, the path that the light must travel is in- 
creased in length, due to the curvature of space-time caused 
by the massive body. This effect is known as Shapiro delay, 
and may be measured for near edge-on binary systems con- 
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J1125-5825 



J1708-3506 



J1731-1847 



J1801-3210 



J1811-2405 








10 cm 




20 cm 



50 cm 




92 cm 



0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 

Pulse Phase 

Figure 2. Average pulse profiles, where available, for each of the new MSPs from observations at wavelengths of 10 cm, 20 cm and 
50 cm using the Parkes radio telescope, and at 92 cm using the Westerbork Synthesis Radio Telescope. In each case, the profile is not flux 
calibrated, and the peak has been placed at phase 0.3. No observations have been made at 50 cm for PSRs J1708— 3506 and J1811— 2405. 



taining an MSP. For low-eccentricity systems, the Shapiro 
delay, Asb, at orbital phase $ is (e.g. Taylor|[l992 1 

AsB = -2rln[l - ssin$] (3) 

where r and s are post-Keplerian orbital parameters, 

r=^Me, (4) 

a = sini (5) 

for an orbit inclined at angle i to the observer with a com- 
panion of mass mc. Defining orbital phase $ = to be at 
epoch of the ascending node, then the Shapiro delay is great- 
est at orbital phase $ = 0.25, when the pulsar is behind the 
companion, and allows both masses and the orbital inclina- 
tion to be constrained where it can be measured. 



( Stairs et al.|2002 l can be calculated, in rads ^ by 



GMr. 



( c3 



V3 



V3 



(mp + TTlc 



73 



(6) 



for a binary with eccentricity e and orbital period Pb. 
The values of calculated for PSRs J1125— 5825 and 
J1708— 3506, assuming the median companion mass given 
in Table [sj are 2 x 10"'' and 6 x 10"^ deg. / yr respec- 
tively. These values compare favourably with th e measured 
cI; = 2.5 X 10"'' deg. /yr of PSR J1903+0327 ( [Champion 
et al.|[2008 l, however, this measurement was obtained with 
a timing RMS of 1.9 ^s, and for a system with e = 0.44, 
which implies that any detection of cj is extremely improb- 
able. 



The magnitude of the Shapiro delay is, however, ex- 
tremely small (sub-/is) for all cases except where the orbital 
inclination is very close to 90° . As shown in Table [s] none 
of the timing residuals for these pulsars are currently below 
5 /IS, so the orbital inclination and the constituent masses 
cannot be constrained. 

Another GR effect that has been observed for some 
pulsars is the advance of periastron, tli. This is best mea- 
sured for orbits with a well-defined eccentricity and angle 
of periastron, cj. Therefore, we only concern ourselves with 
PSRs J1125-5825 and J1708-3506 here. The value of 



3 COMPARISON WITH DISCOVERIES FROM 
PREVIOUS SURVEYS AND PROSPECTS 
FOR MULTI- WAVELENGTH DETECTION 

3.1 Comparison with previous discoveries 

Following the sensitivity curves presented in [Keith et al.| 
( j2010[ ), the left-hand panel of Fig. |3] plots the limiting sen- 
sitivity for the HTRU mid-latitude survey as a function of 
pulse period for pulsars at DMs of 50, 100 and 150 cm"'^pc, 
and also for the PMPS, responsible for the discovery of a 
large fraction of the previously-known MSP population. Also 
plotted are each of the pulsars presented in this paper at 
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Figure 3. Left: Minimum detectable fiux density, for a detection with S/N = 8, for the HTRU survey (solid line) and PMPS (dotted 
line) surveys. For each survey, this limiting flux density is plotted for DMs 50, 100 and 150 cm~^ pc from left to right, assuming the 
scattering law of |Bhat et al.|p004^ . The circles indicate the approximate position of the newly-discovered MSPs in this plane. Right: The 
period and DM of all pulsars with P < 30 ms, excluding those in globular clusters. Squares indicate the position of our new discoveries 
in this plane. 



their period and flux density. Due to the shorter observation 
time of the HTRU mid-latitude survey, the limiting flux den- 
sity at long periods is higher than for the PMPS, however, 
the HTRU survey is more sensitive to pulsars with shorter 
rotational periods, due to the short sampling interval of 64 /xs 
and the high frequency resolution of 0.39 MHz per channel 
(c.f. 250 /is and 3 MHz for the PMPS). The right-hand panel 
of Fig. |3] demonstrates that as well as discovering short- 
period pulsars, the HTRU survey is sensitive to pulsars with 
a high DM. Our new discoveries occupy the short-period and 
high-DM region of the plot, which we can attribute to the 
short sampling time and narrow frequency channels of our 
survey system compared with previous blind searches for 
pulsars (see Keith et al. 2010 for a full discussion). Our 
survey is therefore, on average, finding shorter-period MSPs 
than previous surveys, but are their companions any differ- 
ent? Recently, |BaileS| ( j2007j ) demonstrated that no recycled 
pulsars existed in which the minimum companion mass ex- 
ceeds (P/10 ms) Mq. The only MSP to defy this relation was 
the newly-discovered PSR J1903-I-0327 by [Champion et aT] 
(20081, which is thought to have had a peculiar history. All 



of our pulsars are in agreement with this empirical law. 

In Fig. |4] we plot a histogram of the luminosity at 
1.4 GHz of the previously-known MSPs (where available, 
using P = 50ms as an upper limit), taken from the pulsar 
catalogue ( Manchester et aL]|2005 I . A histogram of the lu- 
minosities of our discoveries, calculated using distances esti- 
mated by the Galactic electron distribution model of |Corde^ 
& Lazio (2002), is also plotted with shaded boxes. We have 
already shown that that the HTRU survey is discovering 
pulsars of lower flux density, at higher DMs. Since DM is an 
approximate measure of the distance to a source, this implies 
that the discoveries have higher luminosities than many of 
the previously known population. This is borne out in our 
results — of the previously-known MSPs with measured lu- 
minosities, ~ 60% have luminosities below 2mJy kpc^. How- 
ever, four of our newly-discovered MSPs have a luminosity 
greater than 2 mjy kpc^ . This trend is expected to continue 



Luminosity at 1 .4 GHz {mJy kpc ) 

Figure 4. Histogram of the 1.4 GHz radio luminosity of MSPs 
(P 50ms) in the pulsar catalogue (plain), including those found 
in globular clusters, and the new discoveries (shaded). 



with future discoveries, providing a more complete luminos- 
ity distribution in the known population. 



3.2 The possibility of detection with Fermi 

In the last year or so, about 30 MSPs have been discovered 
with radio follow-up of non-variable previously unidentified 
point sources from the Fermi Gamma-ray Space Telescope 
( Ray fc Saz Parkinson|20l0 1 . These discoveries form a com- 
plementary group of MSPs to those likely to be found with 
the HTRU survey, since they are faint, but nearby sources. 
A check of the Fermi point source catalogusj^ indicates that 
none of these newly-discovered sources have yet been de- 
tected as gamma-ray sources, after one year of observations 
(it is possible that with further observations Fermi may yet 
make a detection). Compared to pulsars which have been 



http://fermi.gsfc.nasa.gov/ssc/data/access/lat/lyr_catalog/ 



8 S. D. Bates et al. 



discovered from searches of Fermi point sources (Ransom 
et al.|2010 l, this lack of detection is consistent with sources 
which are more distant than the majority of the population. 
A useful metric with which we can compare the likelihood 
of a Fermi detectio n, is by calculating y/E/d? (measured in 



erg kpc s 



Abdo et al. 



2010 



for details). For the 
> 17.5, and 



majority of Fermi-detected MSPs, log(A 
none have been detected where this value is < 16.5. Of our 
discoveries, PSR J1125— 5825 has the highest value of this 

metric, with log(\/i5/ci^) = 16.6. Therefore, it is unlikely 
Fermi will detect this pulsar, or any of our other discoveries, 
assuming the DM distance holds. 



3.3 Potential for optical detection 

In order to decide whether an optical detection of the com- 
panion will be possible for any of the NS-WD systems pre- 
sented here, we must make use of WD evolutionary models 
to estimate the luminosity of the companions. Using WD 



cooling models as described in Fontaine et al. (20011, and 



using the characteristic age of the pulsar as an indication 
of the WD age, we are able to estimate the apparent mag- 
nitude of the companions, assuming the DM distance given 
by |Cordes fc Lazio ( 20C)2| . These models are highly depen- 
dent upon the companion mass, but by using the minimum 
mass (see Table [3|, we are able to put an upper limit on the 
brightness of the companion. Once we include the effects 
of interstellar reddening, it becomes unlikely that the com- 
panion will be visible, each having an apparent magnitudes 
greater than 30. This is a natural consequence of detecting 
sources at large distances in the Galactic plane. 



4 THE ECLIPSING BINARY MSP PSR 
J1731 1847 

PSR J1731 — 1847 is an eclipsing binary system, and only the 
fourth to be published, following B1957-I-20, J2051-0827 
and J1023-I-0038, to be found outside globular clusters. In 
order to understand the origin of the eclipsing material, we 
can approximate the radius of the companion's Roche lobe, 
7?L, as ( |Eggleton|1983[ ) 

Q.maq^''-^ 



0.6(72/3 + ln(l + gi/3) 



0.3 R, 



■0 



(7) 



where q = nic/mp is the mass ratio of the binary system, 
and a is the separation of the pulsar from the companion. 
For this calculation, the companion's mass was assumed to 
be the minimum possible given the mass function, and the 
pulsar mass assumed to be 1.4 M0. In comparison, the ra- 
dius of eclipse is approximately 0.9 Rq, implying that much 
of the eclipse material is not gravitationally bound to the 
companion. 

By studying the delay of pulses near to eclipse ingress 
and egress, we can infer properties of the material causing 
the eclipses. If we assume the delay, at observing frequency 
/, to be purely dispersive, we can express it as a time delay, 
t, relative to a wave of infinite frequency. 



t 



DM 




0.2 0.3 
Orbital Phase 



Figure 5. The bottom panel shows timing residuals for PSR 
J1731— 1847 as a function of orbital phase at observing wave- 
lengths of 10 (-I-), 20 (x) and 50 cm (*). The top panel shows the 
inferred electron density, from equation (|8}, for each TOA. 



Hence, it is possible to attribute a DM to the delay and, 
using the orbital separation of the pulsar and its companion, 
infer the required density of free electrons to cause such a 
delay, assuming a homogeneous plasma. 



4.1 Multi-frequency observations of PSR 
J1731 1847 

In order to probe the eclipse region of this system, obser- 
vations were performed at wavelengths of 10, 20 and 50 cm 
using the Parkes radio telescope. In Fig. [5] one can see that 
no pulses are observed between orbital phases of 0.17 and 
0.35, while the increased residuals either side of the eclipse 
show the extra delay to the pulse arrival times, caused by the 
increased amount of stellar material between ourselves and 
the pulsar as it approaches superior conjunction, which is 
usually defined as phase 0.25. At eclipse egress, one can see 
the additional delay, due to the eclipsing material, reduces 
most rapidly at 10 cm. Hence, using equation (|8|, the added 
electron density can be computed to be 2.5 x 10^^ cm"^ at 
the edge of eclipse. It is also apparent that the ingress and 
egress are not symmetrical about superior conjunction (or- 
bital phase 0.25), with eclipse egress lasting much longer. 
This is indicative of the companion's wind being swept be- 
hind as it moves along its orbit. 

A number of models have been developed to try to ex- 
plain the eclipse phenomenon in pulsar systems, many of 



which are presented in Thompson et al. (19941. The first 



2-KTneC /2 



(8) 



of these is that the eclipse occurs due to refraction of the 
pulse emission by the dense electron content in the compan- 
ion's atmosphere. This model, however, predicts pulses to be 
heavily delayed (typically by between 10 and 100 ms). We 
do not observe this amount of delay, and hence can discount 
the refractive eclipse model for PSR J1731— 1847. Other 
eclipse mechanisms proposed by [Thompson et ah] which 
we disregard in this work are induced Compton scattering 
and Raman scattering of the radio beam. When applied to 
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Table 4. Eclipse properties for the pulsars PSR J1731-1847, PSR J1023+0038, PSR J1748-2446A, PSR B1957+20, and PSR 
J2051-0827. 



Parameter 


J1731-1847 


J1023+0038 


J1748-2446A 


B1957+20 


J2051-0827 


Eclipse Radius (R0) 


0.9 


~ 0.9 


0.8 


0.75 


0.3 


Roche Lobe Radius (Rq) 


0.3 


~ 0.3 


0.15 


0.3 


0.13 


Ph (d) 


0.311134 


0.198 


0.075646 


0.382 


0.0991 


l^'bl 


7 X IQ-" 


2.5 X 10-1" 




3.9 X 10"" 


15.5 X 10-12 


inl/Pb (s-i) 


2.6 X 10^15 


1.5 X 10-1* 




1.2 X 10-15 


1.8 X 10-15 


E/a^ (erglt-s-2 s-l) 


1.5 X 10^3 


3.8 X 10^3 




2.3 X 10^3 


4.8 X 10^2 


Ne Maximum (cm~^) 


3 X IQi* 


4.6 X lO^'' 


2 X 10^** 


4 X lOi'' 


4 X lOi'' 


Anomalous Eclipses 


No 


Yes 


Yes 


No 


No 



PSR B1957+20, these models were found to be unviable 
when explaining eclipses at 1400 MHz, and so we do not 
explore them further. 

Another eclipse mechanism put forward by |Thompson| 
|et al.| is pulse smearing due to the added electron density 
along the line of sight. The first TOA at eclipse egress is 
observed at 10 cm, and is delayed by 450 /xs. By equa- 
tion this implies an added DM of 0.98cm-^pc, which 
would cause the pulse to be smeared by 0.77 ms and 1.3 ms 
at 20 and 50 cm, respectively. This is significantly less than 
the pulse period of PSR J1731 — 1847, and yet the pulses are 
not visible; this indicates that pulse smearing is not wholly 
responsible for the eclipse. 

Absorption of the radio emission by the eclipsing ma- 
terial is another potential mechanism for eclipses, via two 
processes; a) cyclotron absorption; and b) synchrotron ab- 
sorption. However, the optical depth of cyclotron absorption 
is noted by [Thompson et al.| to have a steep spectral index, 
and may not be able to explain the total eclipses seen at 
10 cm (see Fig. |5|. Synchrotron absorption, on the other 
hand, has a shallower spectral index and may be able to ex- 
plain these eclipses we observe in PSR J1731— 1847 at 10 cm. 



these "black widow" systems really are ablating their com- 
panions, this shows some evidence for why this process does 
not occur in all binaries. 

PSRs B1957+20 and J2051-0827 have both been 
demonstrated to show orbital period variations, with both a 
Pb and a Pb. In the case of PSR B1957-I-20, these variations 
take place over a period o f years, with Pb = 1-5 x IQ-n 



and Pb = 1.4 x IQ-i* s-i (Arzoumanian et al. 19941. PSR 



J2051-0827 has a measured value for the Pb of -6.4x 10-^^, 
but this value has been observed to fluctuate over the 14 
years of observations, however, the variations are not smooth 
enough to allow measurement of a P, (Lazaridis, et al., in 
prep.). We do not expect, therefore, to be able to measure a 
Pb for PSR J1731 — 1847 given that our data span less than 
a year. However, we have been able to obtain a tentative Pb 
of -7(2) X 10-11. 

The value of |Pb|/Pb in all these well-studied eclipsing 
systems are of order 10- 1^ 3-1. This gives an estimate of the 
lifetime of these systems of ~ lO'^ yr, although the variation 
of Pb in some of these systems indicates that this lifetime 
may not be too meaningful. From Fig. [l] it can be seen that 
PSR J1731— 1847 has one of the highest P of the eclipsing 
and binary pulsars. 



4.2 Comparison with other eclipsing systems 

Including PSR J1731— 1847 there are now 25 known eclips- 
ing pulsars, many of which are found in globular clusters. 
Table [4] compares the eclipse properties of PSR J1731 — 1847 
with those of four of the most well-studied eclipsing binary 
systems, PSRs J1023-(-0038, B1744-24A, B1957+20 and 
J2051— 0827. All five of these eclipsing systems share similar 
eclipse and orbital parameters, however, PSRs J1023-I-0038 
and J1748— 2446A display anomalous eclipses indicative of 
gas outflow from the companion. PSR J 1023-1-0038 displays 
only a very brief eclipse at ~ 3 GHz, in sharp contrast to 
the case of PSR J1731— 1847. However, this is not unrea- 
sonable since the electron column density, Nc, in the orbital 
system of PSR J 1023-1-0038 (calculated from the additional 
DM given in [Archibald et al.| ([2009), shown in Table |4|, is 



lower by a factor of ~ 6.5. Where Nc ^ 5 x 10 cm- in the 
PSR J1731— 1847 system, the residuals in the lower panel of 
Fig. [5] shown almost no deviation from zero. It is notewor- 
thy that No takes a similar value for each of these systems, 
which is not a requirement of any eclipse models. 

The comparison of E/a^ for these systems shows 
that the typical energy flux at the companion is around 
10^^ erglt-s-^s-i, more than an order of magnitude higher 
than the average for all binary systems. In the case that 



5 DISCUSSION 

The potential for optical detection of counterparts to these 
5 new discoveries is not promising, due to the large dis- 
tance and the associated Galactic extinction. Prospects for 
measurement of Shapiro delay are better, although they 
are dependent upon the orbital inclination. If these systems 
are highly inclined, PSRs J1125-5825 and J1811-2405 wiU 
have Shapiro delays with the largest amplitude, potentially 
with as much as 20 /is delay. In the case of the 6.3 day orbit 
of PSR J1811— 2405, this could be detectable, which would 
allow the measurement of the masses of the pulsar and its 
companion. For all but the highest inclination angles, how- 
ever, the delay will not be observable with the current timing 
precision. 

One of these pulsars, PSR J1731— 1847, is in an eclips- 
ing system with a companion of extremely low mass. Anal- 
ysis of the eclipsing region shows that the pulse delay, if at- 
tributed to dispersion, implies an electron density similar to 
that measured for PSRs B1744-24A and B1957-I-20, which 
are also in eclipsing systems with low-mass companions. 

Two of the pulsars have pulse profiles ideally suited for 
observations with a pulsar timing array. One of these is, 
however, PSR J1731 — 1847, for which there will likely be or- 
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bital period changes on timescales of a year or more, and 
large DM variations that will render it unsuitable for use 
in pulsar timing arrays. As more observations are made of 
PSR J1811— 2405, it will become clear how well precision 
timing can be performed with this pulsar, although the cur- 
rent RMS is not small enough to merit inclusion in a timing 
array. The discovery of distant and bright MSPs is also im- 
portant for a rounded view of the MSP population and their 
distribution in the Galaxy. This has implications for searches 
with future telescopes such as the SKA ( Smits et al.|2009 l, 
and the pulsar timing arrays with these future facilities. 

These discoveries have been made with only ~ 30% of 
the HTRU mid-latitude survey data searched for pulsars, 
indicating that we should expect to find more such objects in 
our data. This compares favourably to simulations presented 
in Keith et al. (20101, which indicated we should expect a 



total haul of ~ 30 new MSPs in the mid-latitude survey. 



6 CONCLUSIONS 

We present the discovery of 5 MSPs in the HTRU sur- 
vey. These MSPs have short periods and DMs among 
the highest in the known population. One of the pulsars, 
PSR J1731— 1847, displays regular eclipses in its 0.3 d orbit 
and with a companion of minimum mass 0.04 M0, appears 
to be a member of the 'black widow' group of MSPs. 

The sensitivity to such objects can be attributed to the 
fast sampling rate and narrow filterbank channels of the 
hardware used for the HTRU survey. 
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